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ABSTRACT 

It is well known that magnetic fields may be represented by the 
cross product of the gradients of two scalars a and ("Euler Poten- 
tials") and that such a representation has considerable use in the de- 
scription of field lines and guiding center motion. Due to practical 
difficulties, however, such potentials are known only for relatively simple 
configurations, including the field of a magnetic dipole. A perturbation 
scheme is developed here by which, given the spherical harmonic ex- 
pansion of the geomagnetic scalar potential and treating it as represent- 
ing a perturbed dipole field, Euler potentials can be derived to any de- 
sired order of accuracy. Its first-order solutions involve trigonometric 
integrals V“ (9) and t“ (6) previously defined by Pennington, which 
may be expressed in closed form; tables of coefficients enabling one to 
obtain these integrals up to n = 6 are included. The results are tested 
for the internally generated geomagnetic field by mapping lines of 
constant a and /3 on the earth's surface; conjugate points (which have 
the same a and ft) are obtained with typical accuracy of 1-2 degrees. 
The method is also applied to a magnetospheric model in which external 
sources are added to represent the effect of the solar wind; although 
this model clearly breaks down at large distances, it gives fairly good 
results within about 10 earth radii, including a "magnetopause" on which 
two neutral points are located. 
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GEOMAGNETIC EULER POTENTIALS 


THE SCALAR GEOMAGNETIC POTENTIAL 

Magnetic fields may be represented by several methods. The most prevalent 
one, and the one most directly connected with the field's sources, is by means of 
a vector potential A. In source-free regions, however, the field may be more 
simply represented as the gradient of a scalar potential, and this is the way in 
which the geomagnetic field is customarily described. We shall denote here the 
geomagnetic scalar potential by y, assuming it to be given by a spherical har- 
monic expansion in spherical coordinates r , 8 and 4> with origin at the earth's 
center (a is the earth's radius): 

B = - Vy (1) 


y - a 


LI 

n=l m=0 



P™ (#) {g™ cos m 4> + h™ sin mcp} 


( 2 ) 


Of the various terms in the expansion of y , the dipole part, corresponding 
to n = 1, is the dominant one. It will be assumed that the polar axis of the co- 
ordinate system has been suitably rotated (giving "tilted dipole" coordinates) so 
that this part is expressed by g° alone, while gj and hj are both zero (Stem, 
1965, Table 1). In this work we investigate quantities which are easily derived 
for the dipole case, but not for arbitrary magnetic fields: it is then useful to re- 
gard the field as that of a perturbed dipole, and carry out a perturbation calcu- 
lation to various orders. 

In any perturbation calculation it is necessary to introduce a "small param- 
eter" e , providing an estimate of accuracy and helping one to decide, at any 
stage, which terms to keep and which to neglect. Since g° exceeds the next 
largest coefficient in the expansion (namely, g* ) by a factor of about 10, the ratio 
between these two terms may be chosen as e . In the perturbation calculation 
one then splits up y into sections of different order 

y = y 0 +y i + r 2 + • • • ( 3 ) 


i 


where y Q is the dipolar part 


7 0 = agj j cos 6 (4) 

and where y contains coefficients ranging between eg® and e 2 g®,y 2 contains 
those between e 2 g® and e 3 g®, and so forth. The perturbation is by no means 
"very small" (as is often the case in perturbation problems), for e is of the 
order of 0.1 and there exist 9 coefficients in the range between e g® and e g®/3. 


EULER POTENTIALS 

The preceding representation suffers from the drawback that it does not 
afford an easy description of magnetic field lines. A representation which allows 
such a description, originally devised by Euler to describe the velocity field of 
incompressible fluids (Euler, 1769, § 26, §49; Truesdell, 1954, § 13), utilizes two 
scalar functions which are conserved along field lines and has the form 


B = Va x 7/3 (5) 

In what follows, such functions will be termed Euler Potentials . They are 
not uniquely defined but may be replaced by any pair u(a, /3) and v(a, /3), pro- 
vided 


3 ( u > v) _ j 

B (a, /3) 


( 6 ) 


For further properties of these potentials, including considerations of their 
single-valuedness and superposition, the reader is referred to more detailed 
work (Stern, 1966). 

Both a and /3 are independent solutions of the linear partial differential 
equation (here written for a ) 


Va • B = 0 


(7) 
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On the other hand, an arbitrary independent pair of functions u and v solving 
this equation does not necessarily satisfy (5). Since the general solution of (7) 
is an arbitrary function of two independent solutions, functional relations then 
exist of the form 

a = a (u, v) 

/ 3 = /8(u,v) 

and substitution into (5) gives 

B = — ■ ~ (Vu x Vv) = )( (u, v) (Vu x V v) (8) 

which, unless (6) happens to hold, is not of the same form. The solutions u and 
v are conserved on field lines and may still be used to label them, but their cross 
product no longer gives the field’s magnitude correctly. Any solution of (7) will 
be called here an Euler potential (the justification of this is given in the next 
paragraph) while a pair of such solutions which also satisfies (5) will be termed 
a matched pair of potentials. 

Given two nonmatching potentials [u, v], it is always possible (at least in 
principle) to derive from them a matching pair [u, S (u , v )] (e.g. Phillips, 1933, 
§49) . Demanding 

B = Vu x VS 

one obtains 

as _ b 

Bv Vu x Vv 

where denominator and enumerator are parallel vectors and therefore have a 
well defined scalar ratio. Consider now the family of surfaces u = constant. On 
each such surface points may be labelled by v and by some other coordinate, 
e.g. 7 or, more generally, the arc length s measured along a field line. Inte- 
grating on such surfaces along lines of constant y (or s ) , from some reference 
point v = v R , one obtains (compare Ray, 1963) 

S = f ? dv' + S'(u) (9) 

J v Vu x V v* 
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Euler potentials have been extensively used in the theory of the motion of charged 
particles trapped by the geomagnetic field. However, there exists no analytic 
representation for them equivalent to equation (2) giving y. The reason for this 
omission is probably the lack of any easy method by which such functions may 
be derived, and it is the purpose of this work to explore means by which this 
situation may be improved, for the case of the geomagnetic field. 


EXPANSION OF THE POTENTIALS 

Because the dipole component of the field possesses axial symmetry, it is 
possible to derive its Euler potentials, namely 

a o = ag ° (f ) sin2 6 (10) 

/ 3 0 = a <6 ( 11 ) 

This suggests that the Euler potentials corresponding to the perturbed dipole 
field can be expanded 

a = a 0 + a-j + a 2 + 

= /3 0 + /3j + /3 2 + 

Substituting these together with the expansion (3) into equation (5) and equating 
individual orders gives 



(Vc^ x V/3 0 ) + (Va 0 x VySj) = -Vy 1 


(13) 


(Va 2 x V/3 0 ) + (Va 1 x V/^) + (Va 0 x V/3 2 ) - -Vy 2 


(14) 


and so forth. Because y 2 represents the limit of accuracy to which the geo- 
magnetic field is determined, there actually exists no practical need for approxi- 
mations of order higher than the second. Each of the above vector equations is 
equivalent to a set of three scalar equations, obtained by forming its scalar 
product with 3 independent vector fields. Let two of these be V a Q and V /3 Q ; the 
equations then obtained are 


(Vaj * Vy 0 ) + (Va 0 • V^) = 0 

(15) 

(VySj • Vy 0 ) + (V/3 0 • Vy 1 ) = 0 

(16) 
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and 


Va 2 

• V y 0 +Va t 

• 

+ Va Q 

• Vy 2 = 0 

(17) 

V/3 2 

• ^7 0 + 



O 

II 

O 

N 

O 

(18) 


These (and corresponding equations of higher orders) are the same as would have 
been obtained from an expansion of (7) to its various orders, and therefore, the 
expanded solutions (12) obtained by considering them alone yield Euler potentials 
which are not necessarily matched. Nevertheless, because the equations in this 
form are linear and involve a 'and $ separately, we shall start by solving them, 
leaving the problem of matching to be considered later. 

It is evident that the functions a t , , a 2 and derived by equations (15) 

to (18) are not unique, for these equations still hold if an arbitrary function of 
a Q and fi 0 is added to any of them. This suggests that we choose a new set of 
independent variables which includes a Q and fi Q . For convenience, we will use 
4> rather than /S Q and retain 6 as the third independent coordinate, even though 
one may more generally consider y 0 or the arc length s 0 along a dipole field 
line for this role. One thus writes 

a = a 0 + a ! ( a 0 0) + a 2 (a Q ,£,<£)+• • (19) 

0 = a4> + (a 0 , 6,4>) + 0 2 (a 0 , 6 , <j>) + * * ( 20 ) 

In these variables, 


d a, i d a, Bar 

Va = _LVa 0 + I -lV/3 0 +— ii- (21) 

1 B a Q u a dcp u do r 

so that 

i By Ba da 3 

V a • V y = _ — B. _ I- = 1 g 0 E-sin Q (22) 

1 ° r 2 ™ ™ a Q<<t> 

where the partial derivative of y 0 , it should be noted, is evaluated with r held 
constant. 

Let Z indicate summation over those terms included in y x . Substitution 
in the equation (15) yields 

Ba \ — 1 n ( dP m 1 

-^-= a/ j < (n + 1) sin 6>P™ (#) + 2 cos Mg™ cos m<£ + h™ sinm<£} (23) 

/in'* ^ 
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Before this can be integrated, the right hand side has to be expressed in terms 
of a Q , 6 and 4>. Using 


one gets 



(24) 



{■ 


sin -2n 6 "S (n + l)sin <9P™ (8) t 2 cos 6 


dP n 


d 9 


(25) 

• {g™ c o s m 0 + h™ s in m 


and hence 


(X\ V” (8) (g™ cosm <b + h™ sin m 4>) + £ ( a 0 . 4>) 


(i) 


agj 


where the function 


(26) 


r e f dp™ i 

V”(0)= sin -2n 0 i (n + 1) sin0P™ (0) + 2 cos 0 j*d© (27) 

*' 77/2 


is a trigonometric integral which may be evaluated in closed form. This result 
is equivalent to one previously derived by Pennington (1966) and tables giving 
V“ (6) are available (Pennington, 1961; Stern, 1965). It also has been derived, 
in a different way, by Birmingham and Northrop (private communication). The 
lower limit of integration has been arbitrarily chosen as rr/2 (dipole equator); 
a different choice merely leads to a different additive function £. 

In a similar manner, the equation for /3 gives 


v— 1 ( a \ n * x 

K 0 l ae 0 J 
( 1 ) \ *1 / 


cos m 6 - gJJ sin m<£} + rj (a Q , <f>) (28) 


with 


r* 

t m (0) = m I sin - ( 2n + x ) 0 P“ (@) d© 

•V 2 


(29) 
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This, too, is equivalent to one of Pennington's results. For the explicit form of 
t™ (<9), see Tables 1 and 2. 


Table 1 

Coefficients of the Trigonometric Integrals t“ ( 6 ) for odd m 
n m 


n+m 

odd 

sin " 1 6 

sin 3 8 

sin- 5 ^ 

sin -7 0 

sin" 9 <9 

sin 1 1 6 



1 

n+m 

even 

Multiply Above Factors by cos Q 

- 

i 

i 

-1. 







2 

i 


-.57735 





.57735 

3 

i 

.73485 

.36742 

-.48990 





3 

3 

-1.58114 

-.79057 






4 

1 



1.10680 

-.45175 



-.65504 

4 

3 



-1.25499 




1.25499 

5 

1 

-.79919 

-.39959 

-.29970 

1.44468 

-.43033 



5 

3 

1.79284 

.89642 

.67231 

-1.79284 

i 




5 

5 

-1.87083 

-.93541 

-.70156 





6 

1 




-2.70044 

2.29128 

-.41660 

.82576 

6 

3 




4.26977 

-2.41522 

i 

-1.85455 

6 

5 

! 

1 



-1.66201 



1.66201 


To derive the second order term a- 2 , it is best to split it into two parts 

a 2 ( a 0 >^> d ) - a 2 + a 2 

where 

(3 a' /3 8 ) gj a 3 r -4 sin 9 - V a Q • Vy 2 (30) 
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Table 2 

Coefficients of the Trigonometric Integrals t™ (#) for even m 


n m 


n+m odd 

sin -2 # sin -4 # sin -6 # sin -8 # sin -10 # 

1 

n+m even 

Multiply Above Factors by cos 

log tan(#/2) 

2 2 

-.86603 

.86603 

3 2 

-.96825 

.96825 

4 2 

.83853 .55902 -1.11803 

-.83853 

4 4 

-1.10926 -.73951 

1.10926 

5 2 

2.56174 -1.28087 

-1.28087 

5 4 

-1.47902 

1.47902 

6 2 

-1.01892 -.67928 -.54343 3.80398 -1.44913 

1.01892 

6 4 

1.39522 .93015 .74412 -2.48039 

-1.39522 

6 6 

-1.25942 -.83961 -.67169 

1.25942 


(3a"/B#) g° a 3 r" 4 sin # = V (31) 

When added together, the above equations yield the equation (17); similar 
expressions may be derived for /3 . 

The equation for a' 2 resembles that for a- 1 and presents no difficulty; in- 
deed, one can easily include the calculation of this part with the calculation of 
a . The equation for a" , on the other hand, requires considerable effort: both 
aj and y x are linear in the set of coefficients (g™ , h” ) summed by I ( x ) and 
therefore the right-hand side of (31) — and consequently a" — contains one term 
for each binary combination of these coefficients. The result can still be inte- 
grated in closed form, but it is rather lengthy and will not be further discussed 
here. 


THE MATCHING OF SOLUTIONS 

The corrections a. and /3. added to the zero order Euler potentials are all 
derived by integration within an arbitrary function of a Q and /3 . The question 
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arises whether the freedom afforded by such a function is sufficient to enable 
one to make the solutions match, within the order to which they are derived. 

It will be shown in what follows — with no attempt at mathematical rigor, how- 
ever — that this indeed is the case. 

Let a subscript such as n characterize quantities of the n-th order, while a 
superscript (n) will refer to a variable or product of variables evaluated to the 
n-th order, i.e. to its perturbation expansion summed up to and including the 
n-th order term. Suppose that matching Euler potentials a (n ) and B (n > of the 
n-th order have been derived by perturbing a 0 and /3 0 , i.e. 


(Va( n ) x V/3( n )}( n ) =_Vy( n ) + 0 (e n + 1 ) 


(32) 


Starting with a< n > and /3 (n) one may then proceed and derive higher orders, and 
we assume that the procedure converges, tending to certain limits a and/3 


a = a.( n ) + a . + a „ + • • • 

n+l T n+2 


/3 = /3< n > +/3 n+1 +/3 n+2 + • • 


In general a and /3, or their approximations to an order higher than the n-th, do 
not match. However, it will be shown that a suitable function k O 0 , fl Q ) can 
be found which, when added to B , causes the (n + 1) order approximation to 
match as well (alternatively, such a function may be derived for a n+J ). The proc- 
ess may then be repeated, allowing matching to any desired order. 

It has been shown in equation (9) that if a and /3 are nonmatching Euler po- 
tentials, a function 8 (a, /3) matching a does exist, i.e. 


V a x V S 


38 

3/3 


(Vax V/3) = - V y 


(33) 


Here we are not as much concerned with 8 as with its derivative 

3 8/3/3 = A. (a, /3) 


We assume that like a and /3, the function A. (a, /3) — and therefore also 8 _ 
can be expanded to various orders 


A. = Xq + + A 2 + • • • 


9 



and so can, therefore, equation (33). Now equation (32) expresses the given fact 
that up to order n , the expansions of a and fi match. When this equation is 
compared to (33), we see that at least to order n, k (a, /3) equals unity 

\("> = 1 + 0(e n+1 ) 


and hence 


^ (n + 1) = 1 +h n+1 (a,/3) 

where h n+ 1 (a,/S ) is some function of order (n + 1). Integration gives, apart from 
an irrelevant function of a alone 


g(n+i) _ ^ + k n+J (a,/3) 

= /3< n+1 > + k n+1 (a , /3) + 0 (e n+2 ) 

Now comes the important step: if one substitutes ( a 0 »/^ 0 ) for (a , /3) in a term of 
order (n + 1) , the error introduced is only of order (n + 2) . In particular 

g(n + l) =£(" + !) + k n + 1) (a 0> /3 0 ) +0(€ n + 2 ) 

This is exactly the desired result, for its states that a solution S (n+1) matching 
a to the (n + 1) order may be obtained by adding to the nonmatching solution 
/3 (n + 1) a function k n+1 (a Q ,j3 0 ) of a Q and . 


MATCHING FIRST ORDER SOLUTIONS 

As has been shown in the preceding section, the first order Euler potentials 
may be matched by a proper choice of the functions € ( a 0 , fi 0 ) and 77 (a 0# ) 

appearing in equations (26) and (28) . It is sufficient for this purpose to consider 
only solutions of the type 




(i) 


ag? 


[V™ ( 6 ) + {g™ cos mcfc + h“ sin mc£} 


/^ = JL (— [t“ (6») + 7J“] {h” cos m4> - g” sin m<£} 

g l ( 1 ) \ ag l / 

where and 17 ” are constants remaining to be determined. 
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The third of the scalar equations equivalent to (13) may be derived by scalar 
multiplication with the unit vector 0 . This gives 


- r 1 ( d 7 j / d #) = ( V / 3 0 x 0 ) • V a.j + (0 x Va Q ) • v p i 


- (a 'r sin 6) (r • Va ( ) - (da Q /dv) (9 • V/^) 


i Ba o J 3a i 
_ + 


r sin <9 3r 1 3a Q 3 0 


or 




a 2 gj sin 3 f9 


oa 3/3 
da o ^ 


(34) 


Expanding, and reconverting the first term to ( a Q , 0, 4>) variables, yields 


£ 


a 0 + * 


ig* con 1.1 -• 4. h“ sin m 4’) sin" 3 0 F® (d) - 0 


n \ 


/ 




dP n 


F m (a) = sin ( 2o ‘ 8 — + n fV m (8) + £ m ] - m [t m (6) + n m l 

n v 7 d 0 n v ' 1 = n L n v/, 'n J 


By making use of Legendre's associated equation it may be shown that 


dF™ 

n 

~de 


•= o 


and therefore F ™(0) is a constant. It only remains to choose the adjustable part 


of F m 

n 


n - m ri m 

^ n 1 n 

so that this constant vanishes. For instance, we may choose 

77™ = 0 
1 n 


li 



and evaluating at 8 = 77/2 we then get 


1 dP; (ir/2) 

n d 0 


(35) 


Table 3 lists the constant parts of [V™ (9) + £“] for n 5 6, with V“ (6) de- 
fined as in (27) and satisfying (35). These values should replace the constant 
parts listed by Stern (1965; Tables 2 and 3) if one wishes to calculate matching 
potentials in the manner described here. 

Table 3 



The above are the ^-independent portions of V™ ( 6 ), selected so that V™ fy/2) = 
(-1/n) [dP™(7r/2)/d 9] 


CONJUGATE POINTS 

In geomagnetic nomenclature, two points (usually chosen on the earth's sur- 
face) are termed conjugate to each other if they are connected by a magnetic field 
line. The identification of such pairs of points is of considerable interest, be- 
cause observations made at them of effects propagating along magnetic field 
lines — e.g. hydromagnetic waves — tend to be correlated (Wescott, 1966). A 
large number of such pairs has been tabulated (Roederer et al., 1966) by numer- 
ically tracing field lines of the geomagnetic field. 

The accuracy of first-order Euler potentials may be tested by deriving a 
map of a and /3 on the earth's surface (Figure 1). On such a map, pairs of con- 
jugate points should be identified by the same a and fi , and the extent to which 
they fail to do so gauges the accuracy of the approximation. 
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LINES OF CONSTANT a AND /3 
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In Table 4 the conjugates to a number of points on the earth’s surface are 
derived using first order Euler potentials (and assuming a spherical earth), and 
they are compared with those derived numerically by Roederer et.al. (1966). A 
deviation averaging about one degree in latitude and a little more than that in 
longitude is evident: it can probably be attributed to the omission of higher- 
order terms. 


Table 4 

Comparison of Conjugate Points 



Latitude 

Longitude 

This Work 

Roederer et.al. 

Conjugate 

Latitude 

Conjugate 

Longitude 

Conjugate 

Latitude 

Conjugate 

Longitude 

Adelaide 

-34.95 

138.53 

52.98 

145.73 

53.2 

143.3 

Anchorage 

61.17 

-149.98 

-53.46 

172.63 

-53.6 

172.2 

Antofagasta 

-23.65 

-70.42 

-4.55 

-68.80 

-3.1 

-69.7 

Archangelsk 

64.58 

40.5 

-52.00 

72.07 

-51.7 

69.1 

Athens 

37.97 

23.72 

-21.25 

31.00 

-19.7 

29.3 

Baghdad 

33.35 

44.38 

-16.53 

48.21 

-16.4 

46.4 

Bangkok 

13.75 

100.55 

5.21 

101.46 

5.4 

100.4 

Belem 

-1.00 

-49.00 

-18.24 

-43.67 

-17.2 

-44.6 

Bismark 

46.82 

-100.77 

-63.18 

133.44 

-62.6 

-135.3 

Bogota 

4.63 

-74.08 

-32.00 

-76.11 

-30.5 

-77.2 

Bombay 

19.00 

72.83 

-0.94 

74.73 

-1.2 

74.3 

Capetown 

-34.15 

18.32 

44.52 

2.83 

45.0 

-0.6 

Deep River 

46.10 

-77.5 

-72.61 

-97.61 

-72.5 

-98.0 

Haifa 

32.83 

35.1 

-15.15 

38.85 

-14.3 

37.0 

Johannesburg 

-26.20 

28.03 

41.38 

19.71 

42.5 

18.3 

Kyoto 

35.02 

135.78 

-18.18 

135.30 

-18.4 

134.7 

Kiruna 

67.83 

20.43 

-58.85 

66.17 

-58.7 

63.6 

Kodaikanal 

10.23 

77.48 

8.41 

78.55 

8.5 

77.6 

London 

51.53 

0.10 

-46.27 

26.14 

-44.3 

25.8 

Midway 

6.97 

158.22 

5.84 

159.10 

6.6 

158.2 

Port Moresby 

-9.20 

147.15 

24.61 

150.36 

24.9 

149.1 

Seattle 

47.75 

-122.42 

-54.18 

-153.73 

-53.7 

-155.2 

Torino 

45.20 

7.65 

-33.67 

23.87 

-32.0 

23.3 

Trinidad 

10.60 

-61.2 

-38.54 

-55.57 

-37.8 

-56.3 

Washington, D.C. 

38.73 

-77.13 

-66.03 

-90.70 

-65.7 

-91.8 
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EXTERNAL SOURCES 


The approximate derivation of Euler potentials by the method outlined here, 
performed in geomagnetic tilted-dipole coordinates (Stem, 1965, Table 1), im- 
proves with increasing distance from the earth, because the field's nondipole 
components fall off faster than the dipole part. However, at distances larger 
than about 3a an additional factor begins perturbing the field, namely the com- 
pression and sweeping of the geomagnetic field by action of the solar wind. It 
is a matter of controversy whether an appreciable current density ("ring current") 
exists in this region; here such currents will be neglected and the external in- 
fluence will be taken into account by adding radially increasing terms to the 
scalar potential expansion 


7 = a 


(i9) {g™ cos m 4) + h™ sin m0} 


n= 1 m=0 


+ a 


n 

LL(ir P™ (<9) {g™ cos m <f> + h™ sin m<£} 


(36) 


n - 1 m=0 


Following Mead (1964), a fair approximation to the external perturbation is ob- 
tained by retaining the terms involving g° and g ‘ along; these coefficients turn 
out to be of the order 10 - 3 gj and 10 -4 respectively. Obviously the field 
can be regarded as a perturbed dipole only up to some (not too well defined) 
limiting distance, of the order of 10a. While the present perturbation method 
will diverge outside this limit, one may use it for the region inside, leading to 
a first order contribution to a 



[V£<0) +£“] 


{g m cos m <$> +h m sin m <b} 

& n n 


where 
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V"(S) 


■ r 

J -n/2 


dp n 


sin 2 ^ n + 1 ) 0 ■< 2 cos 0 — - n sin 0 P m (0)r d0 

1 d 0 nW| 




(38) 



In both cases, the results resemble previously derived ones, with-(n + 1) replac- 
ing n . As before, the constants £ “ and v™ may be chosen so that the potentials 
match. Starting with equation (34), one finds 


a sin" 3 0 



(ag ?/<*„)" {g- 


cos m0 + h™ sin m0} F” ( 9 ) ~ 0 


where 


dP m _ 

F ra ( 6 ) = sin 2n + 3 9 — -- (n + 1) [V m + A m ] - m [t m + 

n v/ A Q v ' n =» n J n =» i 


(39) 


It may again be shown, by means of Legendre’s associated equation, that F” (6) 
is a constant. If if™ andrj™ are chosen so as to satisfy 


(n + 1) 


b n 


+ m "l n = 


(77/2) 


the value of this constant is zero and the matching condition is satisfied. Unlike 
in the case of internal sources, however, additional factors may have to be taken 
into account here when choosing the adjustable constants. 

The reason is that the expansions of a and (3 derived here are valid only 
if Tj and/5j turn out to be relatively small compared to a Q and /3 Q . Because 
of the presence of a Q in the denominators of (37) and (38), this condition will 
not be met, in general, near the poles (9 = 0,n), at which a 0 vanishes. To re- 
duce and possibly remove this source of trouble, it is advantageous to choose 
in such a way (if possible) that [V” + ] vanishes at the poles. From Equa- 

tion (39) it follows that [t™ + 77 “ then also tends to zero, so that such a choice 
also reduces the corresponding divergence in^. 

For m = 0, equation (39) shows that this condition is always satisfied. In 
other cases, may be chosen so that the condition holds at one of the poles, 
e.g. the northern one. ^Whether it then holds for the southern pole as well de- 
pends on the parity of V“(0) with respect to the equator. If (n + m) is odd, 

V”(0) and (9) have both even parity (with appropriate added constant) 

V™(S) =V»(7T-0) 

and hence the condition will hold for both poles. In other cases, the parity (with 
suitable additive constants) is odd and V” ( 6 ) assumes different values at opposite 
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poles (unless it vanishes there). The condition then cannot be met simultaneously 
for both poles: fortunately, this does not occur in the simple model which will 
now be considered, for which only g° and g* are taken into account and all 
parities are even. In that case 

P® = cos 0 P * = /§" sin 6 cos 6 

V° T £ j = - sin 6 6 V* = 2 /3 {(sin 7 6/7) - (sin 9 6/3)} (40) 

t j = Vi = 0 t\ +v\ = (>^377) sin 7 6 

Let G"‘ = g m /gj. The approximate solution is then 

(41) 

a/ag® = -^(a/ r ) - G® (r/a) 2 |- sin 2 8 + 2 73* G * (r/a) 3 {(sin 6/7) - (sin 3 9/3) } cos <fi 
/3/a = cp - (v'3/7) G* (r/a) 4 sin -1 6 sin 4> (42) 


Let A = 10a. Following Mead (1964) we choose 


f® = -0.2511 • 10 -3 
= 0.1242 • 10" 4 


G ® = 0.8 • 10 -3 
G* = 0.3954 * 10 -4 


giving 


(43) 


a/agj = 0. 1 {(A/ r) - 0.4(r/A) 2 } sin 2 6 + 0. 1 (r/A ) 3 (0. 1957 sin/9 - 0.4566 sin 3 6 ) cos </> 


/3/a - cp - 0.09784 (r/A) 4 sin -1 0 sin cp 


(44) 


A plot of a/a g° in the cp - 0,7t plane (corresponding to the noon-midnight 
meridional plane) is given in Figure 2, while Figure 3 gives /3/a in the equatorial 
plane. 


17 



Figure 2— Lines of constant a (which also are field lines) in the noon-midnight meridian of the 
first-order magnetos pheric model. Values given are of \ = ag® /a, which in a dipole field equals 
the equatorial crossing distance, in earth radii, of the corresponding field line. 


It is evident from these figures that up to a distance of about 10a this model 
quite closely approximates our ideas of the magnetosphere's configuration. It 
is therefore not particularly useful to regard it any more, at this stage, as an 
approximation of the 3-parameter scalar potential model from which it was 
derived (and from which it deviates considerably near r = 10a) because that model, 
in its turn, is no more than a simplified representation of the actual field. In- 
stead, we may regard the field represented by (41) and (42) as an independent 
model of the magnetosphere, expressed by means of Euler potentials and con- 
taining two adjustable parameters. The properties of this model — comparison 
to observations, motion of trapped charged particles, etc. — will be explored 
in a separate article. 

We now investigate a little further the meaning of equation (41). Ifg* is 
zero, this solution for a — which will be denoted by a c , the subscript standing 
for "compressed dipole" — is not merely correct to the first order but holds 
exactly, being the well-known solution for a dipole aligned with an external 
homogeneous field (Figure 4). In that case the points at which a = 0 trace the 
z-axis and also cover the sphere 


r = a (2/G°) 1/3 
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: igure 3— Lines of constant /3 (solid) and constant a (broken) in the meridional plane of the 
first-order model magnetosphere. Numbers given are values of A.-ag®/a. 




Figure 4-Lines of constant a (which are also field lines) in the 
meridional plane of a compressed dipole. Letter N indicates 
neutral points. 

which separates field lines attached to the dipole from those reaching to infinity. 
There exist two neutral points where the z-axis meets the sphere, at which the 
field’s intensity vanishes and its direction is therefore not uniquely defined. 

When the first-order contribution of g ’ is added to this configuration, the 
surface of demarcation opens up on one side (the night side) , permitting some 
of the lines originating in the dipole to extend to infinity. The neutral points 
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shift towards the opposite side and the value of a corresponding to them and to 
the demarcation surface now differs from zero, so that the field lines connecting 
them to the dipole are shifted to the front side of the polar field line. The pro- 
perties of the surface of demarcation thus resemble those believed to exist for 
the magnetosphere, the surface observed to separate field lines originating at 
the earth from those originating in the solar wind or in its interaction with the 
geomagnetic field. 

In principle one could perform the preceding perturbation expansion with the 
joint contribution a c of andg° as the zero-order potential. In practice, this 
is not feasible, because the inversion 

r = r (a c , 6) 

which has to be substituted in the expression for 3 a / 3(9 prior to integration can- 
not be readily derived. In addition, the denominator 3 y c / 3<9 of such an expres- 
sion would be too complicated to permit integration in closed form. 
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